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Abstract. Numerical optimization methods such as hillclimbing and simulated 
annealing have been applied to search for highly entangled multi-qubit states. Here the 
genetic algorithm is applied to this optimization problem - to search not only for highly 
entangled states, but also for the corresponding quantum circuits creating these states. 
Simple quantum circuits for maximally (highly) entangled states are discovered for 3, 
4, 5, and 6-qubit systems; and extension of the method to systems with more qubits is 
discussed. Among other results we have found explicit quantum circuits for maximally 
entangled 5 and 6-qubit circuits, with only 8 and 13 quantum gates respectively. One 
significant advantage of our method over previous ones is that it allows very simple 
construction of quantum circuits based on the quantum states found. 
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1. Introduction 

Quantum entanglement [U [2] , enabling states with correlations that have no classical 
analogue, is one of the central concepts differentiating quantum information from 
classical information. Entanglement is essential to many quantum information protocols 
such as quantum key distribution, dense coding, and quantum teleportation, and is also 
thought to play important roles in quantum computational speedup [3]. 

Abstractly, a bipartite pure state is called entangled if it is not decomposable to a 
tensor product of two states of the two subsystems. Quantitatively, there are measures 
[I] which tell us how entangled a state is. Maximally or highly entangled states are of 
particular interest because their entanglement provides a valuable resource which can 
be used to perform tasks that are otherwise difficult or impossible. 

Maximally entangled quantum states of small systems are well known. For 2-qubit 
states, the Bell state |Bell) = ( 1 00) + |11))/a/2 is maximally entangled on all counts: 
maximal entanglement entropy, maximal violation of the Bell inequality, and complete 
mixture of its one-party reduced states. The 3-qubit generalization of the Bell state is 
the GHZ state |GHZ3) = (|000) + \111})/V2, which is also maximally entangled. 

One can easily generalize the GHZ state to general n-qubit states, |GHZn) = 
(|00 - - - 0) =b 1 11 - - - 1))/-\/2. However, for 4 or more qubits, these states are not maximally 
entangled [5J, and in fact display below- average entanglement, as shown by the numerical 
calculations of Borras et al [6] . 

Since the mathematical structures of multi-qubit states are complex, numerical 
optimization methods have been found very helpful in the search for maximally or highly 
entangled states. Using negativity [71 El [9] as the entanglement measure, Brown et al [10] 
performed a numerical search for highly entangled states of 2, 3, 4, and 5 qubits using 
the hillclimbing optimization method. Although it searches through a space including 
mixed states, their method ultimately converges to pure states. They have successfully 
found a simple (in a technical sense of simple to be given below) maximally entangled 
5-qubit pure state, which has application to quantum teleportation, superdense coding, 
and quantum state sharing [TT] . 

There are many entanglement measures besides the negativity [4J. Borras et al [6] 
presented a detailed comparison of the results obtained by different search procedures 
based upon different entanglement measures. They searched only pure states and 
discovered a simple maximally entangled 6-qubit state, the utility of which for quantum 
teleportation and quantum state sharing has been explored [12] . 

Motivated by the simplicity of these maximally entangled 5-qubit and 6-qubit 
states, Tapiador et al [13] conducted a simulated annealing algorithm optimization 
procedure, searching for states not only highly entangled but also algebraically simple. 
They successfully found maximally entangled 5 and 6-qubit states with very simple 
algebraic structure. For 7 and 8-qubit systems, they also discovered highly entangled 
states with rather simple structure. 

The simplicity of the structures of the states mentioned above means [13] that the 
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coefficients of the states with respect to the standard basis of the multi-qubit system are 
nice to write: only a sparse subset of the standard basis vectors have nonzero coefficients, 
and the coefficients are taken from a finite set of simple numbers (±1, ±z, ...). This is 
a natural requirement, since it is harder to derive quantum information protocols based 
on states lacking simple structures. 

However, it is even more natural to demand that the maximally (highly) entangled 
states can be conveniently created by simple quantum circuits. If a maximally entangled 
state can be achieved by a sequence of simple quantum gates, it may be considered to 
have a certain conceptual simplicity, as well as (perhaps more importantly) a superior 
practical realizability. 

Here we report research in which the genetic algorithm is used to search for quantum 
circuits producing maximally entangled multi-qubit pure states numerically. Using 
the negativity entanglement measure, we can construct simple quantum circuits for 
maximally or highly entangled states with up to 6 qubits, a method which, among other 
results, has led us to the simplest known maximally entangled 5 and 6-qubit quantum 
circuits. 

Whether the GA is more effective than other optimization techniques at finding 
maximally entangled states for larger n remains unknown, as the initial computational 
experiments reported here relied upon a GA implementation not suitable for highly 
scalable computation; so exploration of this question remains for future work. However, 
our results do exemplify one significant advantage of our method over previous ones: it 
allows very simple construction of quantum circuits based on the quantum states found. 

Section 2 reviews the entanglement measure used in the paper. Section 3 reviews 
some relevant prior results on maximally entangled quantum systems. Section 4 
describes our use of genetic algorithms, including the details of the encoding scheme 
and the fitness function. Section 5 presents results; and Section 6 gives conclusions and 
discussion. 

2. Measuring Entanglement 

We now review the negativity [3 El [9] entanglement measure which we use as the fitness 
function for our GA. This measure has been shown computationally tractable in several 
previous works pm El [13] . 

The negativity measure is closely related to the PPT [IJJ [15] (positive partial 
transpose) test, which states that positivity of the partial transpose of the density 
matrix of a bipartite state is a necessary condition for the state to be separable. Thus, 
an inseparable (entangled) state is characterized by non-vanishing negativity, where the 
latter is defined as the sum of all the negative eigenvalues of the partial transpose of 
the density matrix. It is convenient to define the entanglement to be the negative 
of the negativity, i.e. = —negativity. The larger E^, the more entangled the state. 
The maximally entangled 2-qubit states are the Bell states having E^ = 0.5. 

It is straightforward to generalize the definition of negativity to multipartite states. 



Constructing quantum circuits for maximally entangled multi-qubit states using the genetic algorithmA 

In this paper, we consider systems with a fixed number of qubits, where each qubit 
constitutes one part of the whole system. Therefore, an n-partite state just means an 
n-qubit state. For an n-qubit system, there are C° + C\ + ... + C™ = 2 n possible cuts 
(partitions). Each cut corresponds to a possible partial transpose. However, since each 
cut has an equivalent cut (e.g. the two cuts {0, 1} [jj and {2,3} for a 4-qubit system 
are equivalent), and we do not need to consider the trivial partial transpose which does 
nothing, there are only 2 n /2 — 1 = 2 ra_1 — 1 nonequivalent cuts. The negativity for an n- 
qubit state is defined to be the sum of the negativity for the 2 n_1 — 1 partial transposes 
and the entanglement Ej$ is still defined to be the negative of the negativity. There 
are upper bounds of E-$ for n-qubit states, which can be derived j6] by considering a 
hypothetical n-qubit state whose marginal density matrices are all completely mixed. 
In this case, each n-cut partial transpose contributes an amount of (2 n — l)/2 to E^. 
For example, the upper bounds of E^ for 3, 4, 5, and 6-qubit systems can be calculated 
to be 1.5, 6.5, 17.5, and 60.5 respectively. See Table 1. 
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3 


4 


5 


6 


1-cuts 


3 


4 


5 


6 


2-cuts 





3 


10 


15 


3-cuts 











10 


total cuts 


3 


7 


15 


31 




0.5 x 3 


0.5 x 4 


0.5 x 5 


0.5 x 6 + 1.5 x 15 




= 1.5 


+1.5 x 3 = 3.5 


+ 1.5 x 10 = 17.5 


+3.5 x 10 = 60.5 



Table 1. Number of cuts for n-qubit systems and the calculation of the (hypothetical) 
maximal entanglement. 



3. Prior Results on Maximally Entangled Multi-qubit Systems 

3.1. Conventions for States, Gates and Circuits 

We use the following conventions for describing quantum states, gates, and circuits. 

The 2" computational basis of the n-qubit states is taken to be (|g n _i)|g n _2) • • • \qo) 
= \ln-1ln-2 ■ ■ ■ %)} where = or 1. A general n-qubit state is expressed as a 
superposition of these basis states as |* n ) = J2q n _ iqn _ 2 ... qo c q „_ iqn _ 2 - qo \q n -iQn-2 • • ■ go), 
where the coefficients c qn _ iqn _ 2 ... qo are complex numbers such that the normalization 
condition J2 qn _ iqn _ 2 ... qo \ c qn -i q „- 2 - qo \ 2 = 1 f° r the state |^ n ) is satisfied. We also denote 
the four Bell basis states as {if)*) = (|00) ± \11))/V2 and |0 ± ) = (|01) ± \10))/y/2. 

The quantum gates acting on a general n-qubit state may be represented by 2" 
dimensional matrices in the above basis. Elementary quantum gates can be conveniently 
expressed as tensor products of small matrices and appropriate identity matrices. For 

§ We label the qubits for an rt-qubit system by the integers 0, 1, n — 1. 
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example, we use the notation H(i) to represent a Hadamard gate acting on the zth qubit 
and use CNOT(i,j) to represent a CNOT gate acting on the zth and the jth qubits, 
which are taken to be the control and the target qubits respectively. 

A quantum circuit of size N will be represented by a string of N elementary gates in 
the form Circuit = Gate(A r -l)Gate(iV-2) • • ■ Gate(l)Gate(0), where Gate(O) acts first 
and Gate(iV — 1) acts last. This convention is chosen to respect the convention of matrix 
multiplication - which is opposite to the convention of the pictorial representations of 
quantum circuits, in which the gate at the leftmost acts first. 

3.2. 4-Qubit States 

Higuchi and Sudbery [5] have proved that there is no 4-qubit pure state with all 
its marginal density matrices completely mixed, which means that the hypothetical 
maximal entanglement (E^ = 6.5) is unreachable. Using variational methods, they 
found a highly entangled state |HS4) with = 6.0981, 

|HS4) = ^(|1100) + |0011)+a;(|1001) + |0110)) + a; 2 (|1010) + |0101)))(l) 
y6 

where to = — 1/2 + iy3/2 is the third root of unity. This state is known to be a local 
maximum [16] and is also conjectured to be a global maximum [5] according to the von 
Neumann entropy measure. 

3.3. 5-qubit States 

Maximally entangled 5-qubit states have been discovered via previous experiments with 
numerical search algorithms [TQj, |6j [13]. These states all have the same entanglement 
feature as the following state, 

|BSSB5) = i(|OOl)|0-) + |010)|Vr> + |1OO)|0 + ) + |111)|^ + )). (2) 

The entanglement distribution of this state among different cuts can be expressed as 
E N = 5 x 0.5 + 10 x 1.5 = 17.5. 

Muralidharan and Panigrahi [11] found various quantum information applications of 
this state and outlined the procedure of a possible physical realization of this state. They 
proposed to create this state by using 2 Hs and 3 CNOTs followed by a 32 dimensional 
matrix with prescribed matrix elements, but the decomposition of this matrix into 
elementary quantum gates was not provided. The advantage of our algorithm, to be 
reviewed below, is that we can find not only maximally entangled 5-qubit states, but 
also the corresponding quantum circuits creating the states. 

3.4- 6-qubit states 

By running a hill climbing algorithm, Borras et al [6J discovered a maximally entangled 
6-qubit state with 32 non- vanishing coefficients, which was used by Choudhury et al [12] 
for various quantum information applications. Maximally entangled 6-qubit states with 
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16 nonzero coefficients were also discovered by Tapiador et al [13] using the simulated 
annealing algorithm. Using our algorithm, we can find both of these kinds of states and 
the corresponding circuits creating them. 

4. Applying the Genetic Algorithm 

We now describe our use of a genetic algorithm to evolve quantum circuits with maximal 
possible entanglement according to the negativity measure. 

4-1. Brief Review of Genetic Algorithms 

The genetic algorithm (GA) [I7J is an optimization algorithm encapsulating the basic 
ideas of biological evolutionary theory (mutation, combinatory sexual reproduction, 
and differential reproduction in a population based on fitness). GAs and other related 
"evolutionary algorithms" have been applied to quantum information and computation 
before, e.g. to the automated generation of quantum algorithms; see [18] and [19] for 
reviews. 

The basic idea of the GA is rather simple. Initially, a starting population 
of individuals (constituting possible solutions to the problem at hand) is generated 
randomly. In our problem, an individual is a quantum circuit which can be used to 
create a new quantum state from an initial state. In the GA broadly conceived, an 
individual has two aspects. On one hand, it has a genotype, which is a sequence of genes 
(also called a chromosome); on the other hand, it has a phenotype, which determines its 
fitness (the quality of the solution). The fitness of an individual is calculated using the 
fitness function which, in our problem is just the entanglement E^- The translations 
from the phenotype to its genotype and vice versa are called encoding and decoding 
respectively. The implementation of these ideas in the context of evolving maximally 
entangled quantum circuits will be detailed in the next subsection. 

After an initial population is created, the GA starts the evolutionary loop, which 
consists of the following steps: 

• evaluate the fitness of each individual in the current population (generation) 

• select individuals based on the fitness levels to be parents of the next generation 

• generate the next generation through genetic operations, namely 

— crossover, a binary operation that takes in two chromosomes and outputs a 
new one combining genes from each of them 

— mutation, which alters a certain percentage of genes in a chromosome, usually 
randomly (sometimes drawn from a particular probability distribution) 

At each generation, the best solution (the elite) in the population is always retained 
for the next generation, until better solutions are found. This loop terminates when 
sufficiently optimized individuals have been found or other predetermined termination 
criteria (such as a maximum number of generations) are met. 
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4-2. Encoding of the quantum circuits 

Our aim, in the present research, is to find simple quantum circuits which can create 
maximally entangled states using the GA. To do this, we first encode all possible 
quantum circuits (solutions) into chromosomes on which the genetic operators can act. 
There are many feasible encoding schemes from which we can choose. Since our problem 
is discrete in nature, we can use binary or integer representations for the chromosomes, 
and we have chosen the latter. 

A quantum circuit is a sequence of elementary quantum gates acting on a number 
of qubits. Any quantum circuit (any unitary matrix) can be approximated by a set 
of universal quantum gates. One set of such universal gates consists of the Hadamard 
gate H, the 7r/8 gate T, and the controlled-NOT gate CNOT. Other useful elementary 
quantum gates include the three Pauli gates X, Y, and Z, the phase gate S, and the 
controlled- Z gate [20]. Not all the gates are needed to construct a quantum circuit 
for a given n-qubit state. For example, the Bell states can be prepared using only a 
Hadamard gate followed by a CNOT gate. Therefore, in our computational experiments 
we have tried various sets of elementary quantum gates to generate states of a given 
number of qubits. 

The elementary quantum gates mentioned above are all 1 and 2-qubit gates. When 
acting on a multi-qubit state, each of these gates can act (nontrivially) on different 
qubit(s). Take the 6-qubit states as an example. If we use H and CNOT as our 
elementary quantum gates, there will be 6 different Hs acting on each one of the 6 
qubits and 30 CNOTs acting on different (ordered) pairs of qubits. We can arrange all 
30 elementary quantum gates into an array, in which each one of the 30 elementary gates 
is specified by an integer between and 35. In this way, a quantum circuit for a 6-qubit 
state using say, 10 elementary gates is represented by an array of 10 integers, each of 
which ranges from to 35. This array of integers is a chromosome, and each integer 
corresponds to a gene. Since there are 10 genes and each gene can take 36 different 
values, the size of the search space (the number of all possible solutions) is 36 10 ~ 10 15 . 
The size of this search space renders brute force search infeasible on current computers; 
but using the GA, we can find an optimized solution using only hundreds to thousands 
of evaluations of the fitness function. 

4-3. Evaluation of the fitness function 

While the chromosome is the object which the genetic operators act on, the fitness is 
the driving force for GA to evolve the population of individuals. The fitness function 
accepts a chromosome as its input and returns a scalar (a real number) as its output. 
The evaluation of the fitness function for our problem takes on the following steps: 

• translate the chromosome (array of integers) into the corresponding quantum circuit 
(sequence of elementary quantum gates) 
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number of gates (3 qubits) number of gates (4 qubits) 




number of gates (5 qubits) number of gates (6 qubits) 



Figure 1. The entanglement generated by the quantum gates at each step of the 
quantum circuits for the 3, 4, 5, and 6-qubit GHZ states. 

• obtain the target state by acting on the initial state (taken as 1 00 • • -0}) with the 
quantum circuit 

• calculate the fitness of the state obtained in the last step using the method 
given above 

5. Results 

We now report the results we have obtained via implementing the above ideas using a 
G A implemented using the genetic algorithm toolbox of MATLAB [2T] , for the case of 
3, 4, 5, and 6-qubit states. 

5.1. 3 Qubits 

For the 3 qubit case, the GA easily finds the GHZ state, as expected. For example, the 
GHZ state |GHZ3) = ( 1 000) + |lll))/\/2 is prepared by the following circuit with only 
3 elementary gates acting on the initial state 1 000) , 



Circuit GH z3 = CNOT(2, 0)CNOT(2, l)H{2). 



(3) 
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number of gates 

Figure 2. The entanglement generated by the quantum circuit for |\E , 4 a (f ) )) at each 
step. 

The entanglement is generated in the following way. Firstly, H(2) change the initial 
state 1 000) to a superposition of 1 000) and 1 100) without generating any entanglement. 
Then, CNOT(2,l) prepares a Bell state (|00) + 1 11)) 10)^/2 for the two qubits q l and 
q2, generating an amount of entanglement E$j = 1. Finally, CNOT(2,0) promotes the 
entanglement to the maximal value = 1.5, making the 3 qubits completely entangled. 
See Figure [U 

Generalizing this construction, the n-qubit GHZ state can be prepared by the 
following circuit with n elementary quantum gates, 

CircuitcHZn = CNOT(n - 1, 0) • • • CNOT(n - 1, n - 2)H{n - 1). (4) 

Again, the Hadamard gate prepares superpositions and the CNOT gates generate 
entanglements. The amount of entanglement generated by CNOT(n — l,m)(0 < 
m < n — 2) is E^ = 2 m_1 and the total entanglement for the n-qubit GHZ state is 
E N = (2™" 1 - l)/2. See Figured] 
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5.2. 4 Qubits 

For the 4 qubit case the result is more interesting. Using H and CNOT as the elementary 
quantum gates, the best solution found by the GA has = 5.5 and the simplest 
quantum circuit we found is of size 5, 

Circuit 4a = CNOT(3, l)CNOT(3, 0)F(3)CNOT(2, 1)H(2). (5) 

The state created by this quantum circuit has the following form, 

|*4a}= i(l 0000 ) + l 011 °) + l 1011 ) + l 1101 )) , fi N 

= ^(|o)l^ + )|o) + |i)|0 + )|i)). w 

There are 7 cuts for a 4-qubit system, 4 single-cuts and 3 two-cuts. The marginal 
density matrices of |^4 a ) for the 4 single-cuts are all completely mixed. For 2 of the 3 
two-index cuts, {0,1} and {0,2}, the marginal density matrices are also completely 
mixed, but the marginal density matrix for the two-index cut {0, 3} is not. The 
entanglement distribution of this state can be expressed as = 4 x 0.5 + 2 x 1.5 + 0.5 = 
5.5. We note that this state has the same entanglement feature as the one used by Wu 
and Zhang to show that not all 4-partite pure states are GHZ reducible [22J . 

By a permutation of the qubits (This is natural, since all the qubits are equivalent 
and there is no privilege for any qubit [23]), this state can be brought to a nicer form, 

l^) = ^(|OO)|^ + ) + |ll)|0 + )). (7) 

The quantum circuit for the state \^4b) can be deduced from Circuity to be 

Circuity = CNOT(3, 0)CNOT(3, 2)F(3)CNOT(l, 0)H(1). (8) 

Although less entangled than the HS state, ^4 a (b) has a simpler form and is much 
more entangled than an average 4-qubit state as can be seen from the entanglement 
distribution of for 4-qubit states [6]. The entanglement generated by Circuit4 (6) is 
shown in Figure EJ Incidently, we note that by using a GA with continuous chromosome 
representation and searching the states only, the algorithm always converges to the same 
entanglement value as the HS state Ej$ = 60.981. This conforms once again that the 
HS state is the global maximally entangled 4-qubit state. 

We can interpret the quantum circuit Circuit4 a (b) in the following way. Firstly, two 
Bell pairs are prepared by two sets of Hadamard and CNOT gates acting on two pairs 
of qubits. The entanglement generated at this stage is E N = 5. Secondly, a CNOT 
gate with the control qubit choosing from one pair of Bell qubits and the target qubit 
choosing from the other pair of Bell qubits increases the entanglement to E^ = 5.5. 

5.3. 5 Qubits 

We have tried various sets of elementary quantum gates in our GA experiments aimed 
at evolution of maximally entangled 5-qubit states. It turns out that it suffices to use 
the set of elementary quantum gates consisting of 5 Hadamard gates H(i) (0 < 4) and 
20 CNOT gates CNOT(z,j) (0 <i,j < 4;i^j). 
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Figure 3. The entanglement generated by the quantum circuit for \^5 a (b)} at each 
step. 



Using these gates as the elementary quantum gates, a maximally entangled state 
similar to |BSSB5) can be easily found by the GA. As we increase the size of the quantum 
circuit, the maximal entanglement obtainable increases monotonically. When the size 
of the circuit equals to 8, the maximal entanglement = 17.5 is achieved. 

One of the maximally entangled states discovered takes the following form, 

l*sa) = Tg(looooo) + loom) + |oion) + |onoo) 

+ |10010) + |10101) - |11001) - |mio)) 

= i(|o)|^ + )|oo) + |o)|0+)|ii) + |i)|r>|oi) + 1 W>|io» 

which is created from the following circuit acting on the initial state 1 00000), 

Circuit 5a = CNOT(3, 2)CNOT(4, l)CNOT(l, 0) 

#(4)CNOT(4, 3)#(4)CNOT(2, 1)#(2). 

By a permutation of the qubits, we can bring the state \^5 a ) to a nicer form, 



(9) 



(10) 



1*56) = -(|ooo>|^+> + |oii)|0 + ) + |ioi)|<r) + |iio)for», 

which can be created from the following circuit acting on the initial state 1 00000), 

Circuity = CNOT(l, 0)CNOT(4, 3)CNOT(3, 2) 

#(4)CNOT(4, 1)#(4)CNOT(0, 3)H(0). 



(12) 
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Figure 4. The entanglement generated by the quantum circuit for at each step. 

The entanglement generated by Circuit 5a (b) is shown in Figure El In this case, two 
Bell pairs are firstly created, followed by a Hadamard gate and 3 CNOT gates. We note 
that at least 3 Hadamard gates are needed to generate maximal entanglement, which 
means that maximally entangled 5-qubit states have at least 8 non- vanishing coefficients. 

5.4- 6 Qubits 

Similarly, we can find a quantum circuit producing maximal entanglement for 6-qubit 
system easily using the GA. 

As with the 5-qubit case, quantum circuits for maximally entangled 6-qubit states 
can be found by the GA using only the Hadamard and the CONT gates. Therefore, 
the set of elementary quantum gates which can generate maximal entanglement for 6- 
qubit states consists of 6 Hadamard gates H(i) (0 < 5) and 30 CNOT gates CNOT(z, j) 
(0<z,j<5;z^j). 

The shortest circuit creating maximal entanglement found by our algorithm has 
13 gates, 5 Hs and 8 CNOTs. A typical maximally entangled 6-qubit state with 32 
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non-vanishing coefficients in this case is, 



^ 6a ) = _^=(|000000) + 1000001) + |000010) - |000011) 

001101) + 1001110) + 1001111) 
010101) - |010110) + |010111) 
011001) + |011010) + |011011) 
100101) + |100110) + |100111) 
101001) - |101010) + |101011) 
110001) - |110010) - |110011) 

liiiiio) + 1 mm)) 



32 

1001100) + 
1010100) + 
|011000) - 
|100100) - 
|101000) + 
1 110000) - 
|111100) - 



(13) 



1 111101) 

= i((|oooo)-|im))(hr) + i0 + )) 

+(|0011) -|noo))(-|v-> + 10+)) 

+ (|0101) + |1010))(|^+) + |0-)) 
+ (|0110) + |1001))(|^ + )- |0-))). 
which is created by the following circuit of size 13 from the initial state |000000), 

Circuity = CNOT(2, l)CNOT(4, l)#(l)CNOT(4, 3) 

#(4)CNOT(5, 2)CNOT(3, 0)CNOT(5, 4) (14) 
#(5)CNOT(3, 2)#(3)CNOT(l, 0)H(1). 

The state |\& 6a ) is maximally entangled, = 6 x 0.5 + 15 x 1.5 + 10 x 3.5 = 60.5. 
See Figure H] for the generation of entanglement of Circuitga- In this case, 3 Bell pairs 
are created before the action of the other 2 Hadamard and 5 CNOT gates. We note that 
\&6a has a more compact form than the maximally entangled 6-qubit state discovered 
by Borras et al 0, [12]. Moreover, by applying a Hadamard gate say, H(0), to \^6a), we 
can get a maximally entangled 6-qubit state with only 16 nonvanishing coefficients, 



(15) 



1*66) = 7g((|0000) - |1111))|^ + ) + (-10011) + |11OO))|0-) 
+ (|0101) + |1010))|V~) + (|0110) + |1OO1))|0+)), 

which is essentially equivalent to the state found by Tapiador et al [13J. 

Our experiments also suggest that 16 is the minimal number of the nonvanishing 
coefficients for the maximally entangled 6-qubit states created from the Hadamard and 
CNOT gates. 



6. Conclusions and Future Work 



In this paper, we have applied the GA to the automated design of quantum circuits for 
maximally entangled states. Highly entangled 4-qubit states and maximally entangled 
5 and 6-qubit states have been found, including the corresponding quantum circuits 
creating them. 

The quantum circuits for the maximally (highly) entangled states are very simple, 
using only 3, 5, 8, and 13 elementary gates for 3, 4, 5, and 6-qubit systems respectively. 
The elementary gates used only include the Hadamard and the CNOT gates. The 
Hadamard gates play the role of generating quantum parallelism (generating non- 
vanishing coefficients of the state) and the CNOT gates play the genuine role of 
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entanglers. Our numerical results suggest that maximally entangled multi-qubit states 
are more-than-two superpositions of the computational basis states. The minimal 
numbers of nonzero coefficients for maximally entangled 5 and 6-qubit states are found 
to be 8 and 16 respectively. This explains partly why the GHZ states are not the 
maximally entangled states for 4 and more qubit systems. 

We note that the construction of quantum circuits for maximally entangled states 
is similar to the construction of perfect entanglers [23] from imperfect ones. A perfect 
2-qubit entangler is a quantum gate which can generate maximal entanglement from an 
unentangled state. The CNOT gate is a perfect entangler for 2-qubit states, but for an 
n-qubit (n > 3) system, the CNOT gates are not perfect entanglers since each one of 
them generates only a (small) part of the whole entanglement of the n-qubit system. In 
this sense, the quantum circuits for maximally entangled states presented in this paper 
can be seen as perfect entanglers for multi-qubit systems. 

The performance of our algorithm appears similar to that of previous approaches 
using hill climbing or simulated annealing. The running time for a single fitness 
evaluation (-En) increases exponentially with the number of qubits. The number 
of evaluations of the fitness function needed before finding an optimal solution also 
increases, though the scaling here is less clear. A typical run of the algorithm for the 
6-qubit case on a common contemporary laptop takes a few minutes, evaluating the 
fitness function several thousand times. 

Envisioned future work includes experimenting with a more scalable C++ GA 
implementation, which may allow search through spaces associated with 7 and more 
qubit systems. We also intend to experiment with alternative entanglement measures 
including the EMM (the entanglement measure based on minors), which is hypothesized 
to be more computationally tractable than eigenvalue based entanglement measures 
|25j . Use of more advanced evolutionary algorithms such as the Bayesian Optimization 
Algorithm (26] may also be helpful; or the use of estimation-of-distribution algorithms 
such as MOSES [27] that act directly on circuit space rather than utilizing an integer 
genome. 
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